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ABSTRACT 

Since the near future should see a rapidly expanding set of probes of the halo masses of 
individual early-type galaxies, we introduce a convenient parameter for characterising 
the halo masses from both observational and theoretical results: V^T, the logarithmic 
radial gradient of the mass-to-light ratio. Using halo density profiles from ACDM 
simulations, we derive predictions for this gradient for various galaxy luminosities 
and star formation efficiencies esp- As a pilot study, we assemble the available V^T 
data from kinematics in early-type galaxies — representing the first unbiassed study of 
halo masses in a wide range of early-type galaxy luminosities — and find a correlation 
between luminosity and V^T, such that the brightest galaxies appear the most dark- 
matter dominated. We find that the gradients in most of the brightest galaxies may 
fit in well with the ACDM predictions, but that there is also a population of fainter 
galaxies whose gradients are so low as to imply an unreasonably high star formation 
efficiency esp > 1- This difficulty is eased if dark haloes are not assumed to have the 
standard ACDM profiles, but lower central concentrations. 

Key words: galaxies: haloes, fundamental parameters, evolution, kinematics and 
dynamics - dark matter 



1 INTRODUCTION 

Spiral galaxies have long been known to show apparent in- 
creases in tlie mass-to-light ratio [M / L) in their outer parts, 
which was one of the prime factors in discovering the ubiq- 
uity of dark matter in the universe ([Bosma fc van der Kruili 
ll979l:IWhitmore. Forbes fc RubiJl988ll . Current studies of 
spiral galaxy rotation curves largely focus on the observed 
radial distribution of the dark matter vis-a-vis the predic- 
tions of the current cosmological paradigm, cold dark matter 
with a cosmological constant (ACDM). In particular, ACDM 
simulations of galaxy halo forma tion predict a steep cusp of 
dark matter in t he halo centres ([Navarro et alj|l996l Il997t 
NFW hereafter; iMoore et all Il999l) . while many observa- 
tions of late- type galaxies indica te relatively low central con- 
ccntrations of dark matter (see'Alam, Bullock & Weinberg" 
[2002: Marchcsini ct al. 2002: do Blok. Bosma. fc McGaueh 
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120031 ICentile et alj |2004l: and references therein). It remains 
to be seen if this discrepancy can be traced to observational 
problems, to oversimplified predictions of the halo properties 
(especially with respect to the inclusion of baryonic effects), 
or to a fundamental failure of the ACDM paradigm. 

Early-type galaxies (ellipticals and lenticulars) are as 
numerous as late-types, and it is important to see how 
their dark matter distributions compare to those of spi- 
rals and to theoretical predictions. Historically, the lack of 
easily interpretable mass tracers like the cold gas in spi- 
rals has made it very difficult to map the mass distribu- 
tion in early-types at large distances from their centres. Ob- 
servations of stellar kinematics from integrated-light spec- 
tra are mostly confined within 2Rc (where the effective ra- 
dius Rc encloses half the projected light), which is gener- 
ally insuf ficient for establishing the presence of a dark mat- 
ter h alo iKronawitter et al]l200o| : iMaeorrian fc Ballantvnj 
I2OOII) — much less for determining its detailed properties. At- 
tempts to constrain the dark matter content in this way 
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demonstrating the difficulties inherent in studying the 
stellar-dominated regions. 

Tracers at larger radius in early-type galaxies are 
needed, and recent advances in instrumentation have dra- 
matically improved the feeisibil ity of measuring halo masses 
using planetary n ebulae (PNe; Douglas e t alJl20o'3) . globu- 

2004), ga seous X-ray emis - 



lar c lusters (GCs; [ Pirsch et^ 



sion jMatsushita et al| 2002': 'O ' Sullivan fc Ponmanll2004bl) . 
gravitational lensing (p<eeton.i200lh . and extended HI disks 
iOosterloo et al ]|2002ir Such studies have so far ascertained 
the presence of both massive ("Romanowsky & KochaM^ 
2001 ; J guote et. al. 2002; Naoolitano ot al. 2002; Goto ct al] 
"2OO3I: iRichtler et alj l2004l : lO'Sullivan k. Pon man 2004b; 
Treu fc Koopmand |2004^ and weak jMendez et aL .2001: 



Romanowskv et al.l200.j|Peng. Ford fc Freemanl2004|) dark 
haloes. Various attempts have been made using halo trac- 
ers to deriv e generic con s traints on dark matter in ellip- 
tical s (Bcrto la et'ai]ll993l:lBahcall. Luhin fc Dornianlll99!Tl: 
IDanzigeB il997l : lLoewenstehi7^^^^Q999ri . but none of 
these has used a large, well-defined sample of ellipti- 
cals covering a broad range of luminosities. With ongo- 
ing observational projects aimed to the study of mass 
distribution around early- type galaxies us ing stellar kine- 
matics (SAURQN: Ide Zeeuw et all l2002t PN.S project: 
3Ugla£_etalJ j2002fl . strong gravitational lensing (CL ASS: 
Myers et al. 1995; and LSD: Koopm ans fc Treull2002ll and 
X-ravs lo'Sullivan fc Ponman ,2004aiilbl) . we expect to have, 
in the near future, a much better observational picture of 
the halo mass distributions in a large representative sample 
of early- type galaxies. 

To prepare for this observational onslaught, in this pa- 
per we want to provide a suite of predictions from the ACDM 
theory for halo mass distributions. To gain qualitative in- 
sight, we begin with the most basic characterisation of a 
dark matter halo around a luminous galaxy: the M/ L radial 
gradient. In conjunction with deriving predictions for this 
gradient, we illustrate the kind of analysis possible with fu- 
ture extensive data sets by comparing the currently available 
observational results with the model expectations. This ap- 
proach may be contrasted with the detailed analysis of indi- 
vidual galaxies, where one attempts to decompose the mass 
distributions into their baryonic and non-baryonic compo- 
nents and then compare these to ACDM predictions. Here 
we are more simply looking at a broad property of the mass 
profiles which can be statistically analysed. Note that even 
with detailed analyses, it may never be possible in individ- 
ual galaxies to unravel all the geometrical and dynamical 
degeneracies, necessitating a statistical approach such as in 
this paper in order to interpret the results. 

The paper is organised as follows: in Section 2 we intro- 
duce the M /L gradient V^T as a basic quantity for describ- 
ing the relative distributions of luminous and dark mass; in 
Section 3 we present predictions from ACDM for V^T, in- 
cluding its dependence on various galaxy properties; in Sec- 
tion 4 we assemble observational results on VfT for a sample 
of galaxies from the literature, and we test for correlations 
with other observed properties; in Section 5 we compare 
the observed and predicted trends in VfT, and examine the 



implications for galaxy formation. We draw conclusions in 
Section 6. 



2 THE M/L GRADIENT 

Aiming for the simplest characterisation of the dark matter 
content of a galaxy, we could consider the total M/L at 
some large radius, such as the benchma rk quantity T b& (B- 
band M/L within 5 Re) introduced bv iRomanowskv et alJ 
l|20o3)- This would allow for constructions similar to the 
Tully-Fisher relations in spiral galaxies, but in comparing 
different galaxies, a quantity such as Tbs is complicated 
by differences in the stellar luminosities (i.e., in the stellar 
M/L, T,). Instead, we consider the increase in the M/L 
with radius relative to its value in the galaxy centre (where 
it is presumably dominated by T*). This M/L gradient gives 
us a clearer idea of how much extra mass resides in the halo, 
and it is computed below. 

The M/L gradient is a continuously varying quantity 
with the (three-dimensional) radius, but its average value 
between an inner and an outer radius (rin,r-out) is given by: 
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where M-t(r) and M^ij-) are respectively the luminous (stel- 
lar) and the dark masses enclosed within a certain radius. 
This is the absolute value of the gradient, but to normalise 
galaxies of different mass and size to similar scales, we com- 
pute 
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Thus VfT is approximately equivalent to the logarith- 
mic gradient of the M/L. The predicted value for V«T is 
the middle expression of Eq. H (also referred to as Vr°'*T 
hereafter) , and is uniquely determined once the mass distri- 
butions and measurement radii are adopted (see Section]^. 
The empirical value is the left expression (V^'^'^T, hereafter); 
in addition to any dynamical modelling uncertainties in AT, 
T» is a source of systematic uncertainty (discussed later). 
Although V^T will in principle depend on the location of 
its measurements (nn/ZJc, rout/i?o), the ratio A/d (r) /A/» (r) 
can, to a first approximation, be assumed to vary linearly 
with the radius (i.e., not to have strong local gradients), such 
that V^T is fairly independent of the measurement radius 
(we che ck this for a the oretical model in Sections l3.3l and l3.4l 
see also lDubinskil998ll . This allows us to uniformly compare 
observational results where the values of {ri^/ Rc, rout / Rc) 
are rather heterogeneous. Note also that V^T is indepen- 
dent of the filter bandpass used for the M/L estimates. 



3 PREDICTIONS FROM ACDM 

We build spherical representations of early-type galaxy mass 
profiles, using a constant- Af/L model for the stellar distri- 
bution (Section 13. H plus a ACDM model of the dark halo 
(Section l3.2ll . We combine these components in Section [3.31 
We do not include a diffuse gas component, since we expect 
its mass inside ^5Re to be a small fraction of the stellar mass 
(we discuss the effect of this assumption in Section I^^J- We 
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derive the predictions for various model parameters in Sec- 
tion OH 

3.1 The luminous component 

The iHernaui'sl il99CJ) profile is known to be a fairly accu- 
rate representation of the stellar mass distribution in early- 
type galaxies; in pr ojection, it approximat es well the de Vau- 
couleurs R^^'^ law Me Vaucouleurjri948t) . Assuming T* to 
be radially constant, the stellar mass density distribution is 

M, k 

P'^"^'' 27V r (r + /fci?e)3' ^ ' 

where k ~ 0.5509. The cumulative mass profile is 

Thus we have a family of density distributions characterised 
by two parameters: the total luminous mass, M*, and the 
effective radius Rc . These parameters do not vary arbitrarily 
in early-type galaxies (as demonstrated by the fundamental 
plane (FP)) but are correlated by a power law of the form 



Rc 



= Rm X 



ft-ilO"A/, 



(5) 



The scatter of the parameters [u^Rm) in Eq. |K| is an 
important source of uncertainty which we will address in 
more detail in later sections. In this section we start drawing 
predictions for generic early-type galax ies in ACDM. For 
this reason we assume the results from IShen et al.l (l200,1 l 
based on a large statistical sample (about 140 000 galaxies) 
as a reasonable representation of the size-mass scaling 
properties of earl y-typ es . 

IShen et alJ (|200.'1l found for their SDSS sample of early-type 
galaxies that a = 0.56 and Rm = 3.6 (Rm = 4.17 for 
h = .7), consistent also with the results of IChiosi et alJ 
l)l998)l . Other studies of large samples of early- type galax- 
ies have found mean relations between total luminosity 
and effective radius of Re oc , with -y = 0.54-0.63 
j|P^rej_^j(^rg2vgld_j^e_(^ gj^ Il998l : iBernardi et alJ 

l2003l: iMamon fc Lokaj|2004a^ . i.e. a slope very similar to 
the IShen et alJ l)200,1) results under the assumption that 
T* is constant with M, (and L). In this case we could take 
a = 7 = 0.6 as an average value. If this assumption is not 
true, as suggested by FP studies, then a relation T, oc Aff 
implies Rc cc L'^ oc Mi^ ''^^ and a — (1-/3)7. The exponent 
/3 depends on the passband and has been foun d to lie in the 
range ~ 0.2-0.6 (Pahrc. Dioreovski & dc C arvalholll995l : 
lOerhard et"ani200ll) : assuming 7 = 0.6 and an intermediate 
value of f} = 0.3, we obtain a = 0.4. Higjher v alues are 
found using the fic ~ Rc relation llKormendvlll977fl and the 
assumption of homology [L oc IcRa, where fic and 7o are the 
effective surface brightness in magnitudes and luminosities, 
respectively). The Kormendy relation ~ 31og7?e implies 
log/e ~ — 1.21og_Rc and with a bit of algebra, 7 — 1.25, 
giving a = 0.9-1.2 for (3 = 0-0.3. 

We thus adopt a — 0.56 as a characteristic value in equation 
5, with a = 0.4-1.2 representing the plausible range. Later 
in the paper, when comparing model predictions with 
observations, we will check that this size-mass relation is 
consistent with our galaxy sample (Section 4.2). Note that, 
besides the relation (a = 0.56, Rm = 3.6) for massive 



galaxies (M, > 2 x 10^°Mq), IShen et al.l J2003l) also 
provide a shallower relation (a = 0.14, Rm = 1.57) for 
dwarf ellipticals. 



3.2 The dark component 

It is a basic prediction of CDM that visible galaxies are 
embedded in extended haloes of dark matter. N-body sim- 
ulations of hierarchical formation in the "concordance" 
ACDM cosmology ( Om = 0.3, Qa = 0.7 ,^8 = 0.9: NFW; 
iBullock et ai]l200ll: IWechsler et~ai]l2002h produce average 
dark halo mass profiles of the form 



Pd(r) = 



(r/r,)(l + r/r,)2' 



(6) 



where rs is the inner characteristic length-scale, correspond- 
ing to the radius where the logarithmic slope of the profile is 
—2. Equation|n]can be written in terms of the total mass 
at the virial radius rvir and of the concentration parameter 



Pd{r) = 
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where Afd ~ 47rps7'f ^(cd) and 

I + X 

The cumulative mass profile is then 



Md(r) = M, 
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where we have fol lo wed the notation of 
iBorriello. Salucci fc Danes3 (|200.tl . In principle, NFW 
profiles are a two-parameter family of density distributions 
(ps and Ts, or equivalently Md and Cd)^. However, a key 
result of ACDM theory is that the halo conce ntration corre- 
latos on average strongly with the virial mass jBuUock et alJ 
^001; Wcchslcr ct al. 2002), which for a population of haloes 
at z — may be written as 
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where the overdensity of the halo mass relati ve to the criti- 
cal d ensity has been taken to be Avir = 101 llBullock et alJ 
l200lfl . In Eq. 1101 Afvir includes the baryonic mass, Mh, 
(Afvir = Md + Mh)- Assuming a baryon fraction /b — 
Qb/fl-M ~ 0.17 (see next section), we have M^ir ~ 1.2Afd 
and Eg. 1101 can be written explicitly for Md as 
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(11) 



^ iNavarro et alj i2004) provide a more detailed analysis of the 
density profiles (see also lMamon fc Loka J2004h|l . but using these 
profiles would change Md out S'Ud V^T by only ~ 10%, which does 
not affect our conclusions in this paper. 

2 We have derived Eq. IIUI in our assumed erg = 0.9 cosmology 
in the halo mass range Mvir = 0.03 — 30 X IO^^Mq via the toy 
model code provided by J. Bullock. 
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Figure 1. Profiles of enclosed mass with radius. Each panel shows a galaxy with parameters indicated. The solid curve is the total 
mass, and the dashed curves are the stellar and dark mass. Insets show the dark-to-luminous mass fraction as a function of the radius. 
These grow almost linearly outside 1 Re- Within 0.5 Re the dark matter fraction ranges between 10-20%. 



Using this correlation, NFW profiles can thus be considered 
as a one-parameter density distribution (characterised by 
Md or Cd). 

The NFW profiles were produced in the frame- 
work of non-interacting (N-body) particles, but real 
galaxy formation includes baryons, whose interactions 
could modify the dark matter mass profiles significantly — 
especially in their centres. Including these baryonic ef- 
fects (cooling, star formation, heating, etc.) can be con- 
sidered the "holy grail" of galaxy formation theory, and 
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|2004). However, it 



seems that even the qualitative effects on the mass profiles 
are not yet known with any certainty, so we are simply tak- 
ing the NFW profiles as the first approximation to the pre- 
dictions of ACDM. Some possible implications of baryonic 
effects are discussed in Section |S] 



3.3 The stellar mass fraction 

From the previous sections, we see that early- type galaxy 
mass distributions in the ACDM framework can be approxi- 



mated fairly well as a family of two parameters (for instance, 
Md and M,). We would further like to parametrise the 
stellar-to-dark mass ratio M*/Md. Although dark matter 
could be stripped from the halo, and baryons could be lost 
or gained through infiows and outflows, for want of more de- 
tailed knowledge we assume that the baryon fraction in the 
galaxy halo is the same as the universal value: /b = fib/f^M- 
This latter quantity is well constrained from the measure- 
ments of the cosmic microwave background, and here we 
adopt /b = 0.17 (Bennett et al. 2003). We relate this to the 
net star formation efficiency esF (including stellar mass loss) 
by the following: 
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and thus esp = 4.9A/, /Afd. While for physical insight, we 
will use esF as our key parameter in this study, it should 
be remembered that this is a valid characteristic only if 
baryons are indeed conserved. Thus we will also quote 
the more robust dark-to-luminous mass fraction parameter 
U=Md/M, =4.9AsF. 

There are no a priori constraints on esF, although 
gala xy formation models use values in the range ~ 0.2- 
1.0 iBenson et alj|20 00'V Its extreme upper limit (given the 
baryon conservation assumption) is esF = 1.0 (/d ~ 4.9), 
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Figure 2. Predictions for the M/L gradient versus star forma- 
tion efficiency. The curves are for stellar masses log(M»/MQ) of 
9.0 to 12.0 (marked in the plot), in increments of 0.5 (stellar lumi- 
nosities Mg of —15.0 to —23.7, assuming T» = 6.5). The top axis 
shows the equivalent dark-to-luminous mass ratio = M^/Mt. 



while its lowest limit can be taken from the universal av- 
erage star formation efficiency o f esp — 0.072 (/d = 68; 
iFukugita. Hoean. fc Peebleslll99^ and the assumption that 
early- type galaxies are at least as efficient as this, since clus- 
ters and groups bring the average down significantly. Stud- 
ies matching meiss and luminosity functions of virialised sys- 
tems have suggested esp ^ 0.25 ( fd ~ 20) for typ ical galaxies 
jMarinoni fc Huds on"200? ; 'Guzik & Scliak'200 j). while stel- 
lar masses in the SDSS survey (Padmanabhan et al.ll20o3) 
implied Mvir/M, = 7 - 30 and /d ~ 6-25 (esF ~ 
0.2 — 0.8). In the following we will consider conservatively 
esF 0.07 — 1.0 as our lower and upper limits respectively. 
We thus construct a sequence of ACDM-based galaxy mod- 
els characterised by M, and esp. In real galaxies, the scatter 
in the R^-M, and Cd-Md relations would produce consid- 
erable variation about this sequence, but for now we are 
considering the mean properties which should be apparent 
in a large enough sample of galaxies. In Fig.0 we show some 
example mass profiles, where typical values for M* are as- 
sumed, and extreme cases of esF = (0.07, 1.0) are used. Note 
that in all cases, the dark matter is the dominant mass com- 
ponent by 5 i?c, indicating that early- type galaxies should 
generically show appreciable M/L gradients by this radius. 
The fraction of the total mass in dark matter within 0.5 Rc 
is in these cases 10-20% (increasing with M* and decreasing 
with esp)'^. It is also evident in Fig.^that the M/L gradient 
is more strongly dependent on esF than on Af* , which arises 
because Rc/ts scales weakly with mass (and thus galaxies 

^ Testing alternative parameter values in the R^-Mt relation 
(Eq. 13 of a=0.4— 1.2 does not change the conclusion that the 
dark matter fraction at is small. 
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Figure 3. Predictions for the M/L gradient versus total stel- 
lar mass from ACDM, for different values of the star formation 
efficiency egp. The solid curves show solutions with the nomi- 
nal Ra-M, relation of (a, Rm) = (0.56, 3.6); they are labelled at 
the right side by the values of egp and the equivalent /d (dark- 
to-luminous mass ratio). Heavy curves are physically plausible 
solutions, and light curves are not. Short-dashed lines show the 
prediction for the shallower R^—Mt relation for dwarf galaxies 
(«, -Rm) = (0.14,1.57): this is only indicative since it is based 
on an i?^/* light distribution, which is not appropriate for dwarf 
systems. Dot-dashed and long-dashed curves show the predictions 
assuming (a, ijjvj) = (0.4,4.4) and (1.2,1.4), respectively; solu- 
tions are plotted for esp = 0.07, 0.25, 1.0. 



with a given /d are a nearly homologous family). In the 
same figure, inset panels show Md(r)/Mt(r); this mass ra- 
tio does not significantly deviate from a linear growth curve 
at r ^ 0.5i?o, except for the case of very low esF (i-c, very 
dark matter dominated haloes) and thus the measurement 
locations in Eq.|21are generally not important (see also next 
section). 

3.4 M/L gradient predictions 

We compute the model M/L gradients V^T (Eq. |5J with 
Tin = 0.5J?c and rout = 4iic representing typical values for 
the empirical measurements in Section r4.1l In Fig.|5| we plot 
the predicted V™°'*T versus esF for sample galaxy models 
with various values of M,. The gradient increases for de- 
creasing esF, since low-efficiency systems imply relatively 
more massive dark haloes. In Fig.|3 we plot the V™°'^T ver- 
sus galaxy mass M* for various values of esp. The gradient 
increases with M, , which can be qualitatively understood 
by the construction that the scale radius of the stars varies 
more rapidly with mass (i?c oc M^'^^) than the scale radius 
of the dark matter (r, oc r^irc'^ oc M^'^^M^ " oc M^ "**^). 
Thus the stellar body of a brighter galaxy encloses a larger 
fraction of its dark matter halo than in a fainter galaxy, and 
appears more dark-matter dominated. This also means that 
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in brighter galaxies, VfT is more sensitive to esp than in 
fainter galaxies (see Fig. I^J. 

Although these predictions appear uniquely determined 
for given (esF,Af«), there is an R^-M, scaling relation im- 
phcitly assumed (Eq. |2J , which translates into assumed lo- 
cations for (rin,rout). We therefore check the effects on the 
model predictions of the Re-M, relation. In Fig.|21we show 
the predicted VfT-M* curves for (Q,i?M) = (0.4,4.4) and 
(1.2,1.4), in addition to the nominal values of (0.56,3.6). 
The steep relation means a steeper trend for V^T: for the 
more massive galaxies, Rc is larger for a given M, , and thus 
the gradient probes farther into the dark matter halo, pro- 
ducing higher values of VfT. For the less massive galaxies, 
7?o is smaller, producing lower values of Vf T. With the shal- 
low relation, Vs scales more quickly than Rc, and the VfT 
behaviour goes in the opposite direction, actually declining 
slightly with M* . 

Note that the mass of any diffuse gas [i.e., Mg = 
M«(l — esp)] should add to the M/L gradient as a bary- 
onic dark matter component. However, assuming that this 
gas is distributed diffusely out to the virial radius (following 
the CDM density profile), its contribution to V^T is no more 
than ~ 0.02. We have also assumed that T» is constant with 
radius, even if observed colour gradients of galaxies suggest 
a slightly decreasing with radius. Although it would be ideal 
to make our measurements in the Jf-band rather than the 
_B-band in order to minimise this effect, we estimate the 
contribution to VfT in our current calculations to be only 
~ —0.02. We also test the sensitivity of V^T to its mea- 
surement locations (rin,r-out); these are important only for 
very dark-matter dominated systems, such that a value of 
V^T ~ 1.0 changes by ~ 0.1 if Arout ~ ±2_Rc. 
Our model galaxy constructions based on ACDM theory 
thus make specific predictions for mass-to-light ratio gra- 
dients. The quantity VfT increases with M, unless esp sys- 
tematically increases with M*, or the scaling parameter a 
is very high. If dwarf ellipticals do indeed follow a shallower 
scaling relation a, then V^ T should increase again at low 
M. (c.f.|Dikir¥sii3ll983), with a minimum occurring at 
logM*//i-iMo ~ 10.2. 



4 OBSERVATIONAL RESULTS 

We describe our compilation of AI/L results for early-type 
galaxies in Section 14.11 We study the Rc-M, relation in 
Section f4.2l and search for correlations between parameters 
in Section r4.3l Hereafter we assume Ho — 70 km s~^ Mpc^^. 

4.1 Galaxy sample 

We assemble from the literature the dynamical results on 
all early-type galaxies with mass measurements at a radius 
''out i5 2_Rc. These exclude galaxies with masses derived 
through X-ray analyses, strong gravitational lensing and gas 
disks and rings, as these techniques may be prone to strong 
selection effects. For instance. X-ray emission occurs mostly 
in giant elliptical galaxies, so these results are restricted to 
the brightest, most massive s ystems — mostly group or clus- 
ter dominant ellipticals (see lO' Sullivan fc PonmanI l2004bl 
for a discussion). Gravitational lensing also suffers a selec- 
tion effect for large masses, since the lensing cross-section 



is a power-law of the lens mass. Furthermore, the systems 
studied by l ensing are generally loca ted at large redshifts 
0.1 < z < 1 iTreu fc Koopman jliooi) : here it appears that 
low mass galaxies are still formin g stars while only very mas- 
sive galaxies are already in place jjimenez et al I2OO4) . This 
makes the selection effect even worse if we want to com- 
pare results with local galaxies. Extended HI rings (which 
are rare occurrences) are mostly located in SO systems. Fur- 
thermore, in most of the cases, these systems are gas rich 
{Mm/ Lb.q ^ 0.2) which suggests that these galaxies are 
the result of a particul ar formation history ( merging rem- 
nants of gas rich disks, lOosterloo et al 1 1200211 and are not 
representative of the ellipticals on the whole. Thus, none of 
these methods is really suitable for an unbiased survey of the 
properties of dark matter haloes around elliptical galaxies. 

The techniques remaining include integrated stellar 
light, PNe and GCs; there is no obvious reason for any of 
these to produce a particular selection effect on total mass 
(rather than on luminosity). The results could in principle 
be skewed by dust e ffects when relying on i ntegrated stel- 
lar light kinematics iBaes fc Deionghgl200 j) . but the halo 
masses of most of the sample galaxies are derived from dis- 
crete tracers (PNe and GCs), which are not affected by this 
issue. There are 21 galaxies in the sample (16 ellipticals, 3 
lenticulars and 2 of a transitional nature); their properties 
are summarised in Table Q including 04 (isophote shape pa- 
rameter) and 7 (inner surface brightness slope). Note that 
this sample includes two cD galaxies, whose halo masses may 
very well be dominated by the galaxy cluster in which they 
reside. However, this may also be true of group-dominant 
galaxies in general, and at this point we do not wish to 
make such a priori distinctions; we will check any density 
environment dependence later on. 

Some of the dynamical mass-to-light ratios (Tin, Tout — 
provided in the B-band for uniformity) are derived 
from highly flexible modelling techniques which fit or- 
bital distri butions to the detailed li n e-of-sight velocity dis- 
tributions llKronawitter et al.l I2OO0I: iGebhardt et al.1 120031 : 
lRomanowsk^^^ni2003ll . For these we take the uncertain- 
ties in T as given in the papers. In other cases, Jeans (and 
similar) models have been used to fit binned velocity disper- 
sion profiles; in such models, orbital anisotropy is a major 
source of systematic uncertainty. Although not many de- 
tailed dynamical studies have extended beyond 2 Rc, so 
far n o cases of extreme a nisotropy have been detected (see 
e.g., iGerhard et aIll2001^ . We therefore assume that the 
anisotropy parameter /3a lies in the range ±0.3 (/3a = is 
isotropic), and estimate uncertainties in Tin of ±10% and 
in Tout of ±30%. 

To find the empirical gradient V^'^'T, we also need to 
know T,, which may be estimated in two completely inde- 
pendent ways. We can use spectroscopic information for the 
luminous body combined with (single-burst) stellar popula- 
tion synthesis models: T* = Tssp. This is in principle the 
most direct way to measure the stellar mass only, but cur- 
rently such models are susceptible to considerable system- 
atic uncertainties. The second way is to use the dynamical 
estimate: T* = Tin; there will be some dark matter contri- 
bution to Tin, so this approach will provide a lower limit on 
yobsy rpj^g advantage of this approach is that many sys- 
tematic uncertainties (such as distance) will cancel out in 
Eq.|5| the empirical gradient is thus: 
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Table 1. Catalogue of early- type galaxies with halo mass measurements. 



Galaxy 


Typc*^ 


(Mpc) 






^4 


d 


1 








^out 
(Rc) 




^SSP 




Data type ; 
'^in,out I^^f- 


NGC 221 


E3 


0.9 


-16.5 


0.2 












0.1 


2.8 ± 0.2 


5.6 


4.8 ± 2.8 


„ „-|-0 9 


0.13 ± 0.18 


LS,PN;MB01,V02 


NGC 821 


E2 


25.5 


-20.6 


6.2 


2 


5 





64 


0.5 


8.4 ± 0.4 


4.8 


13.1 ± 3.9 


+ 1.9 

— 2 


0.13 ± 0.10 


LS + PN;G03.R03 


NGC 1316 


SO 


22.7 


-22.1 


12.0 












0.6 


3.2 ± 0.3 


2. 


5.7 ± 1.7 


3 3 + 0. 8 


0.56 ± 0.39 


LS + PN;A98 


In K-y 1 -J 


EO 


21 2 


19 8 


2 5 










4 


4 4 ± 4 


2 ' 


4 4 ± 4 


7 4+2-^ 


zb 07 


L S ' M B 1 


NGC 1399-^ 


El/cD 


21.1 


-21.4 


4.3 





1 





07 


1, 


8.3 ± 0.8 


12 


42 ± 13 


10 4+'^ '^ 
-2.7 
3.3 ± 0.4 


0.37 ± 0.14 


LS + GC;S00,R04 


NGC 1700 
NGC 2434 
NGC 3115 


E4 
El 

E6/S0 


52 
22.8 
10.2 


-21.7 
-20.5 
-20.2 


3.5 
2.7 
2.8 




1 


4 









01 
7 

78 


0.5 
0.5 
0.8 


4. ± 0.4 
8.2 ± 0.3 
6 ± 0.6 


4.6 
2.6 
3.3 


7.8 ± 0.8 
15.8 ± 3.4 
10 ± 2 


0.23 ± 0.06 
0.44 ± 0.20 
0.27 ± 0.14 


LS;S99 
LS;K00 
LS;C93,E99 


NGC 3379 


El 


11.2 


-20.1 


2.0 





2 





18 


0.5 


5.9 ± 0.4 


7.9 


8.7 ± 1.1 


-■-til 


0.06 ± 0.03 


LS + PN;G03,R03 


NGC 3384 


E5/S0 


12.3 


-19.8 


1.5 


1 








6 


0.5 


4.1 ± 0.4 


5.8 


7.4 ± 2.2 


3.2 ± 0.6 


0.15 ±0.1 


LS,PN;T95,B96 


NGC 4406 


E3 


18.1 


-21.6 


8.9 


-0 


7 





08 


0.3 


6.4 ± 0.7 


2.4 


11.4 ± 3.4 


-1^21 


0.37 ± 0.23 


LS,PN;M91.A96 


NGC 4464 


so/a 


15.3 


-18.7 


0.45 










88 


0.5 


7.9 ± 0.8 


2.8 


8.7 ± 2.6 


-IS-.; 


0.04 ± 0.15 


LS;MB01 


NGC 4472 


E2 


17.2 


-22. 


8.7 


-0 


3 





04 


0.5 


8. ± 0.3 


4.5 


28.5 ± 8.6 




0.64 ± 0.27 


LS,GC;K00,C03 


NGC 4486 
NGC 4486B 
NGC 4494 


E3/cD 

El 

El 


17. 
17. 
IS 


-21.7 
-16.8 
-20.7 


7.8 
0.26 
4.3 







8 
3 







2 

14 
6 


0.5 
0.2 
0.5 


5.3 ± 0.4 
8.1 ± 0.6 
3.9 ± 0.4 


4.8 
2.8 
3.9 


30. ± 4.5 
10.3 ± 2.2 
5.5 ±1.7 




1.1 ± 0.2 
0.1 ± 0.11 
0.12 ± 0.13 


LS + GC;RK01 
LS;KOO 

LS,PN;K00,R03 


NGC 4697 


E3 


12.4 


-20.2 


5.7 


1 


4 





74 


0.4 


9.3 ± 0.4 


3.4 


9.3 ± 2.8 


6.7+-3 


0.0 ±0.1 


LS,PN;A'I01 


NGC 5128 


SO 


4.4 


-21. 


6.6 


-0 


5 





15 


0.3 


2.8 ± 0.3 


10 


10 ± 3 




0.27 ± 0.07 


PN;P04 


NGC 5846 


El 


26.3 


-20.4 


8.0 







<0.2 


0.5 


9. ± 1. 


3.2 


20 ± 6 


4.7 ±1.0 


0.45 ± 0.25 


LS;MB01 


NGC 6703 


El 


28.2 


-20.3 


4.1 











0.5 


5.3 ± 0.3 


2.6 


7.6 ± 2.3 


0.21 ± 0.21 


LS;KOO 


NGC 7626 


E2 


46 


-21.4 


8.5 





2 


<0.2 


0.5 


10. ± 0.6 


2.1 


18.3 ± 3.6 




0.52 ± 0.23 


LS;KOO 



Types and magnitudes from HyperLeda fhttp: / / www-obs.univ-lyonl .fr/hypercat/1 . 
" Distances from SBF measurements <Tonrv et all200j) . rescaled from h = 0.74 to h = 0.70. The exceptions 
are NGC 4486B— which is assigned the same distance as NGC 4486— and NGC 1700, NGC 4464 and NGC 
7626, for which the redshift-distance is used. The distance uncertainties are taken from the SBF paper and 
included in any uncertainties in Mg and Af, , but are not included in the T uncertainties shown here since 
they cancel out when calculating V^T. 

*^ For consistency. Re is taken from the dynamical reference paper. 

a4 an d 7 are ta ken fro m |Benderj.^urstcin fc Faber' |'19 92^: iForbes . Franx. fc lUingworthl ^^9^; 
iFaber et al. 11997); Pcl legriniUl99Hh : lTrager ct al. 12000); Ravin dranath et a l. ( 2001); and references therein. 
^ LS=long-slit integrated light stellar kinematics; PN=planetary nebulae; GC=globular clusters 
/ The PN and GC resuhs from .Napolitano et al.i Ii200a) {M/L ~ 28-45 at 8.8 Rc) imply V^T^ 0.3-0.58 
depending on assumptions about the halo equilibrium. We use the R04 result when computing correlations. 
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We first check the reliability of these two approaches 
by comparing their results in a common galaxy sub- 
samp le;_TlM_jnean_stel^ age and metallicity are taken 
from [ T erlevich fc Forbe'3 ( |200^, or if unavailable , from 
iTrager et alJ J2000l) oTfrom iKuntschner et'aP (1200 ll) : note 
Tin = ^0/8 generally. We then obtain the expected B- 
band stellar mass-to-light ratio, Tssp, usi ng stellar popula- 
tion synthesis models from lWorthevI il994) . After examining 
galaxies in common in order to tie the results onto a uniform 
scale, we increase the Tssp val ues from Trager ct al. (2000i') 
by 15% and decrease those from IKuntschner et al.l ll200ll) by 
8%. The results for the 15 galaxies with available age and 
metallicity are listed in Table Q We compare Tssp to the 
innermost dynamical estimates Tin in Fig. |1] and find them 
to be remarkably consistent given all the potential sources of 
error which have not been included in the uncertainties (e.g., 
distance, IMF, heterogeneous Vin/Ra)- This consistency in- 
dicates that if some fraction of Tin comes from dark matter, 
then Tssp has been systematically overestimated by approx- 
imately the same amount. We will thus assume T, = Tin 
for all our ensuing analyses, also examining in due course 
the potential effects of a non-negligible dark matter contri- 
bution. 

A best fit of T, to galaxy luminosity and mass shows 



a weak correlation: Tb.* oc L^o^+o.oe^ ^ ^ ^^o.ioio.oe^ 
This is weaker than some of the trends found in th e lit- 
erature (sec Section | 3. Ill .but iKauffmann et aP ll2003l) and 
^admanabhan ct al. ( 2 0o3) also find T« to be roughly con- 
stant. The models of the latter two imply T, 5 in the 
B-band, while we obtain ~ 7 typically, lending credence to 
the possibility that we are underestimating the amount of 
dark matter in the central parts of our galaxy sample. 



4.2 Mass-radius relation 

As seen in Section 2, to derive the M/L gradient predictions, 
we needed a mean relation between Re and M* (Eq. 01, 
where we adopted values for the parameters (a,i?A/) based 
on literature results. For direct comparison to our empiri- 
cal AI/L gradient results, we must check that our galaxy 
sample is consistent with this relation. We begin by fitting 
a mean R^-Mt relation to our sample (as shown in Fig. 
where M* is obtained by multiplying T* and the total lu- 
minosity Lb- We find {a,RM) = (0.50 ± 0.13 3.6 ± 0.8) 
which is consistent with (0.56, 4.2; h = 0.7) from lShen et alJ 
1.2003') . Since this fit is sensitive to sparse sampling effects 
in our data set, we instead wish to use an exponent a based 
on the much larger SDSS galaxy sample, with the normal- 
isation Rm set by our data because of systematic differ- 
ences in measuring masses and sizes. Assuming a typical 
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Figure 6. Empirical M/L gradients versus galaxy parameters: the total luminosity M-q (upper left), the total stellar mass Af» (upper 
right), the isophotal shape 04 (lower left), and the central surface brightness slope 7 (lower right). Error bars are l-cr level. 



a — 0.6 (see Section we find Rm = 3 .1 ± 0.8; for 

comparison Rm = 4.2 from Ish en et alJ (|200^ and the 
IPahre. Diorgovski fc de CarvalhoTl^g^ results (where also 
a ~ 0.6) imply Rm = 2.6 if {B - K) = 4.1 and T. = 7. 
Thus, our Ra-Mt relation is intermediate between typical 
literature values. The uncertainties in this relation are im- 
portant to consider when making predictions for V^T. In 



* This is consistent with their T* values (obtained via population 
synthesis modelling; Kauffmann et al. 2003^ being systematically 
lower than ours. 



particular, the exponent a is not yet known conclusively 
and may systematically vary with M, ; we have seen in Sec- 
tion l3.4l how changes in a would affect the model predictions 
for VfT. We have also seen that, even if a is known accu- 
rately, there is a large scatter in real galaxies around the 
mean R^-Mt relation (which is after all only a projection 
of the fundamental plane). So when we later come to M/L 
results in Fig. |7| these should not be over-interpreted for 
any single galaxy; rather, the mean trends in a large sam- 
ple of galaxies are meaningful. We can quantify the effects 
of the scatter by finding an "envelope" for the Re-M, rela- 
tion, enclosing 68% of the data points. The upper and lower 
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Table 2. Spearman correlation test 
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Figure 4. Estimates of stellar M/L from stellar population syn- 
thesis modelling (Tssp) compared to dynamical modelling (Tin), 
in units of Tg^Q. The dashed line shows Tin = Tggp. 
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Figure 5. Relation between total luminous mass and effective 
radius. The solid line is the best fit to the data (o = 0.6, Rm = 
3.2) and the dotted and dashed lines show the ± 1 cr uncertainties 
of the fit: (a = 0.74, _RAf = 4.0) and (a = 0.46, ^Af = 2.4), 
respectively. 
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N 


^^spcar 


confidence level 


V^T - Mb 


21 


-0.58 


99.9% 


VfT - 04 


21 


-0.42 


98% 


ViT-7 


19 


-0.39 


97% 


V^T - M, 


21 


0.58 


99.9% 


V^T - Pcnv 


18 


0.10 


78% 



envelopes have {a,RM) = (0.74,4.0) and (0.46,2.4), respec- 
tively (see Fig. |KJ ; these values are well within the plausible 
bounds on a of 0.4-1.2 discussed in Section f3.1l Thus when 
we compare to theoretical models in Section 5, we can safely 
use (a = 0.6, Rm = 3.1), based on larger galaxy surveys, as 
representative of our sample; and we can check the effects 
of allowing (q,_Ra/) to vary within plausible limits. 



4.3 Correlations 

Before comparing the empirical results on the M/L gra- 
dient to theoretical predictions, we would like to see if 
there are any correlations in the data (see earlier work in 
.Capaccioh et al. 2003; NaDolitano ct al. 200i). The broad- 
est property to test against is the total galaxy magnitude 
Mb- Plotting V?'""T against Mb in Fig. ^ we see that 
fainter galaxies have shallow gradients (V^'^'^T ~ 0.0-0.1 for 
Mb ^ —20) while the brighter galaxies have a wide range of 
gradients (Vf^'T ~ 0.1-0.6). The same trend is found when 
considering galaxy mass M« (as estimated using T» — Tin) 
rather than galaxy luminosity, with logM* ~ 11.2 marking 
the transition between shallow and steep gradients. We also 
see evidence in Fig. |S] for disky galaxies (04 ^ 0.2) having 
shallow gradients while boxy ones have a wide range, and 
for galaxies with a steep central cusp (7 J; 0.2) having shal- 
lower gradients than those with a shallow cusp. 
These last two trends are not surprising given the first, since 
it is well-established that galaxy propertie s such as a4 and 
7 correlate strongly w i th luminosity (Nieto fc Benderlll989l : 
ICaoaccioh et alJll992t iFaber et al^ll997^ . The strengths of 
the co rrelations as meas ured by the Spearman rank statistic 
' spear llPress et alJl992^ are given in Table 2, where A'^ is the 
number of galaxies. To estimate the statistical significance 
of these correlations, we performed, for each structural pa- 
rameter, 10 000 random experiments where V^T have been 
randomly extracted for each galaxy according to a Gaussian 
measurement distribution. The confidence level of the cor- 
relations being real is then given in the table. 
The correlations apparent by eye are all significant at the 
more than 97% level, and the ones with Mb and M^, are the 
strongest. Note that these correlation results are indepen- 
dent of any distance uncertainties (which cancel out in the 
method used for estimating V^'^'T). Nor are they changed if 
Tt is systematically wrong in the sense that there is a large, 
constant fraction of dark matter in the galaxy centres: this 
would change V^'^'^T equally for all galaxies and would not 
affect their ordering. A possible effect is if T» = Tin is wrong 
by an amount which varies systematically with Mb- 

However, in our basic ACDM picture, T« is likely to be 
systematically too high for the brightest galaxies (see Sec- 
tion leading to V°'''^T estimates which are too low, 
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Figure 7. M/L gradients versus total luminous mass. Points 
with error bars show galaxy data with uncertainties. Labelled 
curves are our predictions from ACDM for different values of esF 
and /tj: solid curves are models for the best-fit i?c— M* relation 
(a = 0.6, Rm = 3.6), short dashed and long-dashed curves are 
models for the 68% envelopes as in Fig. ISl 
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Figure 8. Efficiency of star formation versus total luminous mass, 
for early- type galaxy data fitted to ACDM models (open squares). 
Plausible upper and lower bounds on these parameters are shown 
as dotted lines. Full points are the median values from the distri- 
bution of the best-fit values. The data points of NGC 2434 and 
NGC 4486 are above the plot limits. 



thus, the "corrected" values of °T would result in an 
even stronger correlation. Even if T* were wrong in the op- 
posite sense (such that T,/Tin is actually much lower for 
low-L ellipticals), we estimate that dark matter mass frac- 
tions of ~ 50% would be needed to explain the correlation — 
much larger than the fractions typically estimated from 
other studies. 

As a final source of systematic effects in the M/L trends, we 
have checked any dependence of V^'^'^T with the (density) 
environment. We used for this purpose the galaxy d ensity 
estima tes penv (number of galaxies per Mpc'') from ITuUvI 
l|l988l) and compared with our V^^'^T. We did not find any 
significant correlation with M/L gradients: the Spearman 
rank statistics gives a 78% significance for a positive corre- 
lation. The steepest gradients (NGC 4486 and NGC 4472) 
are indeed the ones related to the densest environment (but 
note that their close companions NGC 4486B and NGC 4464 
have low gradients) . These are the only apparent cases where 
the environmental effect could have enhanced the gradients 
because of a cluster dark matter backgrou nd or the effect 
of a close encounter (see iNapolitano et aiJ l2002l for NGC 
1399). 

We therefore conclude that there is a empirical trend for 
bright, boxy shallow-core early-type galaxies to appear more 
dark-matter dominated than faint, disky steep-core galaxies. 
This result is irrespective of the type of underlying dark mat- 
ter distribution and environment, but is especially strong if 
it resides in basic ACDM haloes. 



5 IMPLICATIONS OF PREDICTIONS VS. 
OBSERVATIONS 

In Sect ion [5. II we show the results on halo masses and im- 
plied esF using the V^'^'^T data from Section0]in conjunction 
with the model predictions of Section f3. 41 Finding many of 
the halo masses to be problematically low in the context of 
ACDM, we consider possible alternative explanations in Sec- 
tion In Section [5.31 we briefly discuss the results based 
on other observation techniques (X-rays, gravitational lens- 
ing and HI rings) and we lastly give some indication of the 
esF inferred in a few systems from the very low-mass regime 
in Section [5.41 



5.1 Matching models and observations 

We now compare the galaxy M/L gradient data, V^'^'^T, to 
the expectations from ACDM, and illustrate this in Fig. |71 
where Vf^'^T is plotted versus M* for various values of esp 
and using the best-fit Rc-M^, relation (a = 0.6, Rm ~ 3.1) 
together with the 68% envelopes. We can see immediately 
that the data do not seem to follow the simple increase of 
VfT with M^, that is expected for a "universal" esp, and 
the empirical increase appears to be sharper. By inspection 
of the figure, one interpretation is that the low-mass galax- 
ies have very high values for esp (in some cases even above 
the highest allowed limit) while the high-mass galaxies have 
a broad scatter in esp, mostly within the allowed efficiency 
range. Note that projection effects in the dynamical mod- 
elling of T can produce systematic effects on the inferred 
yobsy. fa^e-on disky galaxies could appear to have very low 
Tout if interpreted as spherical systems. This would produce 
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Figure 9. M/L gradients versus total luminous mass, with 
curves shown of constant concentration parameter cj, according 
to the ACDM predictions, with egp free to vary. Line styles are 
as in Fig. |7| 



some data points with spuriously low V^'^'^T (and high esp) 
values, but these should comprise a small fraction of a large, 
unbiased sample. 

To quantify the esp-Af* trend requires allowing for the 
individual galaxies' departures from the mean Rc-M, rela- 
tion. To this end, we fit the model simultaneously to Rc and 
to V° "T with free parameters (a, Rm, £sf), minimising the 
following statistic: 



(14) 



where qi = Fi = i?M(M.)", g2=Vf=T, 

F2 =V™°'*T(a, _Rm, esp) and Sqt are the uncertainties on qt. 
The results are shown in Fig.|H]where we plot the best-fit esp 
(and equivalently /d) against the stellar masses. From the 
distribution of the fitted values we have inferred the follow- 
ing median values: esp ~ 2.0ln a for logM* — ll.O^n'^ and 
esP ~ 0.23t° " 



-0.8 ^"s - 
-0 13 ^o'' ^''S ~ ll-4!to 2 (with a marginal over- 



lap of the mass ranges around log M* ~ 11.2). There are also 
galaxies with logM, = H-Tto l for which esp ~ 0.041°;°^, 
only marginally consistent with the plausible lowest effi- 
ciency limit. These values suggest a smooth relation between 
M* and esp (see Fig. |HJ which is at this stage only tentative 
because of the large uncertainties in the data. 
An alternative interpretation is the presence of a mass scale 
[log(Mo/M0) ~ 11.2, Mo ~ 1.6 X lO^M©] marking a clear 
dichotomy: galaxies with M, ^ Mo are most consistent with 
esp ~ 0.25 (/d ~ 19.5), while galaxies with M* ^ Mo are 
most consistent with esp ~2 (/d ~ 3). In the latter sam- 
ple, half of these are well determined to have esp > 0.9 
(/d < 5.4), while those galaxies with lower esp (larger V^'^'^T 
in Fig. |7J are boxy systems, in agreement with the tentative 
dichotomy discussed at the end of Section 14.31 This sug- 
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Figure 10. Baryon fraction for the low-mass systems (Af* < 
Mo) as a function of the halo concentration Cd: here the high 
Cd ^ 13 range is shown. Halo masses are inferred by their best fit 
Cd according to the Cd — and shown on the upper axis. Dashed 
line marks the cosmological baryon fraction, /jj = 0.17. 



gests that there are separate populations of galaxies with 
very different esp, rather than a trend of esp with M« for 
a single galaxy population — but we cannot distinguish be- 
tween these two scenarios given the large uncertainties in 
the current data. 

Note that other studies have suggested a maximum 
galaxy formation efficiency (minimum dar k-to-luminous 
mass fraction, /d) at M« close to our Mn iBenson et all 
L2OO.0I; iMarinoni fc HudsonI I200I iRomano et alJ 12002^. ~ 
suming an average T« = 6.5, as in our sample, we can con- 
vert this mass scale into a luminosity scale of Mb = —20.5, 
fairly similar to the luminosity scale found in other lumi- 
nous properties of early-t ype galaxies iGraham et alj|2003l : 
iGraham fc Guzmanlliool) . 

The esp ~ 0.25 value for the high-mass galaxies is loga- 
rithmically midway between the uppermost and lowermost 
limits adopted. However, esp ~ 2 is unphysical for the low- 
mass galaxies, given the assumptions of our ACDM models. 
To understand the origin of these esp values in the ACDM 
models, in Fig. El we plot VT°'^T again as a function of M* 
with curves of constant Cd now indicated. At each point on 
the curve, esp has been adjusted to fit the required ACDM 
Cd-Afd relation. For a low stellar mass galaxy (Af» < AIo) 
to have a shallow AI/L gradient, its ACDM concentration 
must be high (cd 5; 15), which implies a small halo mass. 
This mass corresponds to esp > 1 and thus to a violation 
of baryon conservation. To show this, we plot in Fig. 1101 
the baryon fraction /b as a function of Cd, where the /d re- 
sults (Fig. (HJ have been re-interpreted with the assumption 
esp = 1.0, giving /b — 1/(1 -I- 1.0/d). We find that many of 
the galaxies with Cd ^ 15 have /b ^ 0.17, a circumstance 
which get worse if esp < 1. 
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Figure 11. M/L gradients versus total luminous mass account- 
ing for 30% dark mass in the inner regions. Line styles have the 
same meaning as in Fig. 13 The position of the data points rela- 
tive to the model curves is almost unchanged with respect to Fig. 
Ill 



Figure 12. M/L gradients versus total luminous mass and con- 
centration parameter after c^—M^ relation relaxed for esF = 0.25 
. Labelled solid lines are the predictions from the ACDM for dif- 
ferent C(j: line styles as in Fig. 171 



5.2 Discussion 

In the previous section we have seen that the basic assump- 
tion of ACDM (namely the Cd — Md relation) and baryon 
conservation with /b = 0.17 lead to typical formation ef- 
ficiencies esF ~ 0.25 for galaxies with M, > A'Iq, while in 
the low-mass regime (M* < Mo) they produce unphysical 
esF • We have also concluded that if we want to decrease the 
inferred esF below LO, and still have /b = 0.17, we must re- 
lease the Cd — Md relation otherwise we produce a violation 
of the baryon conservation. 

Before moving on to other possible explanations for this 
problem, we first examine whether this problem can be ex- 
plained by the systematic uncertainties in T,. As discussed 
in Section 13.31 in the ACDM framework, we expect ~10- 
20% of Tin to be a dark matter contribution. Reducing T« 
by as much as 30% will systematically shift all the empirical 
estimates of V^'^'T higher by 30% (see Eq.Ell. The R^-M* 
relation will also be affected, such that Rc for a given M* 
is higher, and thus rin/rj, and rout/rs in the model calcula- 
tions will be higher; this leads to a higher value of V™°'*T 
predicted by the models which roughly offsets the higher 
values in the data. The final effect is for the data values 
of M* to shift, and so the net change is for the low-mass 
galaxy data points in Fig. |7|to shift roughly parallel to the 
contours of constant esF — not affecting the conclusions on 
esF significantly (see Fig. IIH . To decrease the inferred esF 
for the intermediate-mass galaxies to values below 1.0 would 
require a dark mass fraction inside Tie of ^ 50%. Such high 
central dark matter fractions are not supported by other 
empirical studies of the central mass content of early-type 
galaxies llCerhard et al.lf200ll : iRusin. Kochanek. fc KeetonI 



1200311 . They would also imply that the average stellar mass- 
to- light ratios, Ts., = 6.5 ± 2.0, must be decreased by a 
factor of two: this is not plausible from the standpoint of 
stellar populations modelling where T_b,* ^ 4 only for fairly 
young stellar populat ions, with age ^ 5 Gyr, and met allicity 
-0.5 < [Fe/H] < 0.5 iWorthevlll99i : lMarastonlll998^ . 
One way to avoid the problem of low-mass haloes is to relax 
the Cd-Md relation, using a maximum plausible value of esF 
as a prior. Once esF and M* are fixed, the total dark-halo 
mass, Md, is set and the only way to alter Md(r) (Eq.jnj is 
by altering Cd. In this way V™°'*T can be written as a func- 
tion of Mt and Cd: in Fig. 1121 this is shown for esF ~ 0.25 
(/d = 19.5). We can derive the best-fit Cd values for each 
galaxy using the same procedure adopted for best-fitting 
esF (as in previous section). 

Results are in Fig. 1131 (left): we see that the most massive 
galaxies (Md ~ 8 x IO^'^Mq) have concentrations consistent 
with ACDM expectations (cd ~ 8-14), while many of the 
less massive galaxies (Md 2 x 10^^ Mq) have very low con- 
centrations (cd ~ 1-6). 

The shallow Ad/L gradients in the low- mass galaxies are thus 
a simple consequence of a small fraction of the dark matter 
re siding in the galaxy centre. This explanation was advanced 
in lRomanowskv et alj i2003l) but contradicts ACDM expec- 
tations that Cd is decreasing with Md rather than increas- 
ing. As suggested by Section [5.11 esF may not be constant 
but rather decreasing with galaxy mass, which may allevi- 
ate the problem. We thus assume esF=0.25 for M* > Mo 
and esF = 1-0 for M* < Mo (a maximum plausible value), 
and plot the results for Cd in Fig. 1131 (right). The Cd-Md 
relation is now in a better agreement with the ACDM ex- 
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pectation, but the low-mass galaxies still have haloes with 
concentrations which are too low on average. 

The above alternatives (baryon non-conservation; low 
concentrations) for producing low AI/L gradients are ne- 
cessitated by the strict NFW functional form Md (r) for the 
haloes (Eq.|UJ. Thus there is little flexibility to change Md,out 
vs Md,in- But another way to reduce V^T in the models 
(Eq. |5J is to change the ratio of Md,out versus Md,in by 
changing the form of the mass profile Md(r). For example, 
in the classic adiabatic contraction model of galaxy forma- 
tion, the baryons' collapse drags dark matter inward; this 
will increase Md,in above the nominal NFW halo central 
mass and decrease the inferred Vf^^T. Again, the implica- 
tions of this alternative would be that many other studies 
of central DM density are incorrect, and that the popula- 
tion synthesis models have over-estimated T, by ~50%. As 
mentioned in Section f3. 21 it is not clear that this adiabatic 
contraction scenario is really to be expected, and we leave 
examination of the details to a future study — but comparing 
the alternatives, it seems this may be the most likely one for 
"saving" ACDM. 

Of course, the ACDM theory is not as readily 
falsifiable as may be supposed, as various modifica- 
tions to its simplest outline had been proposed in or- 
der to solve various problems. One may assume a dif- 
ferent initial power spectrum and in principle partially 
solve the problems seen here; it has also been suggested 
that (Tg ~ 0.75 is a reas onable solution to the low- 
concentration halo problem l|M cGaug h. Barker fc de Blokl 
l2003l:lvan den Bosch. Mo fc YandbOOaV 

5.3 Comparison with other dynamical techniques 

For the sake of completeness, in this paragraph we briefly 
discuss how the results based on stellar kinematics plus 
PNe and GCs compare with results listed in the literature 
based on the techniques we have excluded in this paper: 
namely X-rays, gravitational lensing and HI rings. This 
is to verify whether (regardless any possible selection 
effects) different techniques give consistent results for 
the dependence of the M/L gradients on the stellar 
mass and luminosity, which we have explained as a trend of 
the star formation efficiency with the same global quantities. 

X-rays. lO'SuUivan fc Ponma 3 ^2004b^ reviewed litera- 
ture results on M/L estimates based on X-ray emissions 
in early- type galaxies. They found that at 5-6 Re different 
modeling techniques give a lower limit of M / L=2Q Mq /Lq 
on a large galaxy sample. Even if most of these estimates 
come from cluster and group dominant galaxies, they 
concluded that ellipticals can possess dark haloes of their 
own, regardless of their environment. We can infer a lower 
limit for the V^'^'^T for these systems using their results. 
We can assume a Tin = 7 as the average value find in our 
sample and derive V°''''T(Xrays)> 0.4. Since Mb ^ -20.5 
in the X-rays sample, this is consistent with esF ^ 0.25 
found for our galaxy sample of stellar masses > Mq. 

Strong Gravitational lensing (SGL). There are only a 
few studies where M/L of the innermost and outer regions 
are available for reliable V^'^^T estimates. As discussed 
above, there is the further complication that SGL analyses 



refer to high z {z ^ 0.1) and they cannot be directly 
compared to model predictions obtained in the local uni- 
verse {z = 0) since both dark and luminous matter scaling 
relations change with redshift. A detailed model prediction 
at high-z is beyond the purpose of this paper. Here we want 
to g ive only a qualitative inter pretation of SGL results. 

iTreu fc KoopmansI (l2004ll reviewed 5 bright {Mb ^ 
—20.5) E/SO systems with reliable M/L estimates inside 
the Einstein radius (ranging from 0.8 to 5 Re)- Using 
their results we compute V^'^'^T^ 0.4 and an average 
V?''°T(SGL) ~ 0.8. These values are consistent with but 
somewhat higher than our dynamical results for bright 
galaxies, and comparison to z = models implies a very 
low star formation efficiency (esF ^ 0.1). 

iRusin. Kochanek. fc KeetonI i2003l) used statistical 
constraints on the mass profiles of 22 lens systems over a 
broad range of luminosities to find a typical homologous 
mass profile. Their results assuming a ACDM dark matter 
profile imply V^T ~ 0.1-0.2, which is consistent with our 
res ults for faint galaxie s but a pparently not consistent with 
the iTreu fc' KoopmansI (1200411 results and with our bright 
galaxy results. A more careful examination of these results, 
including treatment of the z-dependence of the models, 
awaits a future paper. 

HI rings. For these galaxies, the gradient prediction 
should account the presence of the gas since its mass 
contribution it is not negligible for these systems. For the 
sake of simplicity, as a first approximation we use our 
predictions obtained without g reference. There are 

a few systems in the literature with reliable innermost and 
outer M/L values suitable for V°'^'^ T estimates. 
These g alaxies include NGC 5266 llMorganti et al.lll997ll . 
IC 5063 fMorganti. Oosterloo. fc T svetanov"l99^. IC 2006 
(|Pranx, van Gorko m. fc de Zeeuwl.l 994V NGC 1052, NGC 
2974 and NGC 4278 iBertola et alJll993t) . These galaxies 
span a broad range of luminosity, and we derive Vf =T 
values ranging from 0.4 to 1.0, corresponding to esp ^ 0.2. 
These results are consistent with what we have found 
for the brighter galaxies, but the weak halo behaviour 
for fainter galaxies is not detected. However, as discussed 
earlier, we consider it dangerous to take the HI results as 
typical of early-type galaxies because of selection issues. 

5.4 Dwarf galaxies 

In Section|3we have made V^T predictions for a wide range 
of galaxy masses, including dwarf systems (logAf* <^ 10). 
This mass regime is populated both by dwarf elliptical 
galaxies (dEs) and low surface brightness late-type systems 
(LSBs). These two classes are different in many respects 
(morphology, gas content, stellar population etc.), but they 
also share some similarities: they are small, low-luminosity 
galaxies with diffuse, exponential declining surface bright- 
ness pro files^. Furthermore they have the same size-mass 
relation iShen et al. 2003?) and hence they are expected to 
have similar dark-to-luminous mass distributions. Specifi- 
cally, as shown in Section we must expect observed gra- 

^ In this respect our basic model assumption of an R}^'^ surface 
brightness profile is poorly representative for these systems. 
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Figure 13. Best fit values after relaxing the cj-Md relation for esF = 0.25 (left) and esp varying within the galaxy sample (right: 
see details in the text). Solid line is — from BOl with l-cr scatter indicated by dashed lines. 



dients to increase for decreasing stellar masses (and lumi- 
nosities) and for decreasing star formation efficiency (at any 
given mass) . Unfortunately, these galaxies remain among the 
most poorly studied systems, mostly because of the observa- 
tional difficulties related to their faint luminosities and de- 
tailed dynamical analysis are scarce and generally focussed 
on th e central galaxy regions jBender fc Nietdll990l : iMated 

Eiii). 

A complete view of these systems, although beyond the pur- 
pose of this paper, is complicated by data which are gener- 
ally not suitable for the kind of analysis based on M/L gradi- 
ents. Here we want to briefly compare our models with a few 
literature data sets which can be adapted to this analysis. 
We restrict our analysis to dwarf early-type galaxies which 
are gas poor and can be directly compared to the predictions 
of Section|Hl where we discarded any gas component. 

Fornax JWalcher et al.l2003^ [Ms =-11. 7]: they estimate 
M/L=23 at 90' from the galaxy centre {Rc = 9.9') and a 
central M/L^A.S. We have obtained 'T=0.52 for this 
galaxy which corresponds to esF=0.6 using the dwarf shal- 
lower Re — M« relation in our models. Wo found consistent 
esF values using M/L estimates from Lokas (2001). 

FS373 llDe Riicke et al.1 12004H [Ms=-16.9]: they ob- 
tained M/L=7.8 at 1.5 7ie(=7.9"/1.57 kpc) and computed 
a stellar M/L=2-A from age and metallicity of the stellar 
population. This allowed us to estimate V°''''T~l-2.5 and 
esF=0.1-0.6; 

FS76 iPe Riicke et alj |2004 [Ms=-16.1]: they ob- 
tained M/L=7A at 1.5 i?e(=4.4"/0.77 kpc) and considered 
a stellar M/L=2-4 as for FS373. We found pretty similar 
V?'''T~l-2.5 but esF=0. 2-0.45 for this galaxy. 

This very preliminary test on dwarf ellipticals suggests 
that M/L gradients are actually increasing with decreasing 



masses as expected, and that this trend is also consistent 
with esF tending to decrease at this mass scale from a max- 
imum near L* . 

Low star formation efficiencies in this mass range 
have been also inferr e d for gas-rich LSB galaxies 
Jjobii^^^angnan[ Il99d: ICote Carignan. fc Sancisil Il99ll : 
iMartimbeau. Carignan. fc Rovl Il994^ . Thev found L ~ 
Mtot/M* =25-27^ in the galaxies NGC 3109, NGC 5855 and 
IC 2574 respectively, corresponding to observed esF ~ 0.2 
which are almost in agre ement with efficiencies of dwarf 
ellipticals derived above. iLake. Schommer. fc van GorkornI 
|l99(|) claim at least half of the baryonic mass of DDO 
170 in the form of gas, which implies esF ^ 0.5 assuming 
baryon conservation. All these esp estimates are consistent 
with what we obtain by comparing the V^^'^T from their 
M/L estimates with our models (in spite of our model sim- 
plification of not including gas). 

We can then conclude that, based on AI/L gradients, there 
is evidence of low star formation efficiencies in the very low 
mass regime. This s eems in agreement with o ther studies 
iBenson et alJl200(]l : iMarinoni fc HudsonI l20o3) claiming a 
continuous decreasing star forming efficiency with decreas- 
ing masses at these mass scales. 



6 SUMMARY AND CONCLUSIONS 

We have introduced a method for measuring AI/L gradients 
in early-type (elliptical and lenticular) galaxies, quantified 
by V^T. If the same technique is used for all M/L measure- 
ments in a galaxy, then VfT is independent of bandpass 

^ Mium in their paper refers to the total baryonic matter 

-^-^(stars-t-gas) . 
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definitions and distance uncertainties. We have built typical 
galaxy models in the ACDM framework using NFW profiles 
for the dark matter distribution. Assuming baryon conserva- 
tion across the virial radius, these galaxies are parametrised 
by their total stellar mass M, and their net star forma- 
tion efficiency esr; we have then made predictions for V^T 
(between 0.5 Re and 4 Re) as a function of (M*, esp). As a 
consequence of the relative spatial scaling relations of the lu- 
minous and dark matter, we have found that brighter galax- 
ies should typically show higher values of V^T than fainter 
ones, appearing more dark-matter dominated even though 
the overall dark-to-luminous mass fraction /d is a constant. 
These predictions are valid for a mean set of galaxies only, 
and should not be used to draw direct conclusions about 
individual galaxies, whose properties are subject to consid- 
erable scatter. 

We have next assembled a set of AI/L measurements 
(using stellar, PN and GC kinematics) in 21 early-type 
galaxies as a pilot study for comparing to our model predic- 
tions. We have found that V^T in these galaxies increases 
with M, (and with other correlated parameters), which is 
evidence for systematic changes in either the fraction of dark 
matter or in its radial distribution. This empirical trend 
must be explained by any successful model of galaxy prop- 
erties and formation (including those advocating alternative 
theories of gravity) . Environmental density does not appear 
to play a key role in producing these trends. 

Comparing the empirical V^T values to the ACDM pre- 
dictions, we find a stellar mass scale Mo ~ 1.6x lO^^M©, and 
a corresponding luminosity scale of Afs ~ —20.5, marking 
a galaxy dichotomy. The brighter (more massive) galaxies 
{Mb ^ —20.5) show a broad scatter in esF, but are gen- 
erally consistent with esp ~ 0.25. The fainter (low-mass) 
galaxies are generally only consistent with esp ~ 2-3, which 
is not physically possible in the context of our model as- 
sumptions. 

The brighter galaxies thus appear consistent with the 
ACDM picture without need for further complications from 
baryonic effects. For the fainter galaxies, we have examined 
some possibilities for explaining their apparent extremely 
high esp values (i.e., apparent extremely low dark matter 
fractions). One explanation is that the baryon fraction has 
not been conserved in these galaxies: either they have ac- 
quired large quantities of extra baryons from infiows across 
the virial radius (e.g., gas from nearby starbursting galax- 
ies, which has since formed new stars); and/or they have 
lost portions of their outer dark haloes (after star forma- 
tion) through dynamical interactions with other galaxies and 
groups. We note that the latter process is implicitly included 
in ACDM simulations of dark halo formation, but the net 
effect on /d has not been explicitly studied. 

Another explanation is that the predicted relation 
between halo mass Md and concentration Cd in ACDM 
is wrong, and Cd tends to increase rather than de- 
crease with Md, with Cd ~ 1-6 for L* galaxies. Simi- 
lar conclusions have already been reached in studies of 
low-surface brightness galaxies ("Alam. Bu llock fc Weinberd 
[2002 : McCa ugh. Barker & de Blok 2003), and may indicate 
that the dark matter is not CDM, or that the baryonic pro- 
cesses involved in galaxy formation have somehow caused 
the dark haloes to become less centrally-concentrated. How- 
ever, the current best guess for baryonic effects is that 



the cooling and collapse of the gas cau ses the dark mat- 
ter halo to also contrac t adiabatically jBlumenthal et all 
119861: iGnedin et ahlliool . which increases the central den- 
sity of the halo. This is still a viable scenario for the fainter 
early-type galaxies if the stellar AI/L values are systemat- 
ically lower by 50% than what we have assumed, and the 
central dark mass fraction is ^40%. This possibility merits 
more investigation, although we note that several studies 
have concluded that the dark matter cannot be so domi- 
nant. Both scenarios (low-concentration haloes; low stellar 
M/L) are possible ways for resolving an apparent empiri- 
cal inconsistency: weak lensing and satellite dynamics stud- 
ies find evidence for large amounts of dark matter inside 
the virial radius in L* ea rly-types JCuzik fc Seliakll2002l: 
Ivan den Bosch et alJl2004^ ■ while this paper suggests little 
dark matter inside ~ 4f?c in similar galaxies. 

Whatever the case, these data imply a dichotomy in 
the dark halo profiles between fainter and brighter early- 
type galaxies (see Section 14.31 . which parallels the di- 
chotomy alre ady visible in the properties of their lumi- 
nous bodies (Nieto fc Bonder '1989"; CaDaccioli et al. 1993; 
Fabcr ct al. 1997; Graham ot a L 2003; G r aham fc Guzi naj] 
2003), and X-ray properties llPellegrinil Il999l) . and sug- 
gests markedly different formation histories for the two 
galaxy populations where baryon cooling efficiencies could 
have played a pivotal role. The comparison of data and 
models also supports other studies indicat ing that esp 
is at a maximum for galaxi e s around L* iBenson et 
hood iGuzik fc SeliakI 120021: iMarinoni fc HudsonI l2002l : 
Ivan den Bosch. Mo fc Yang|l2003^ . Indeed, we have briefly 
discussed the case of dwarf/LSB galaxies and shown that 
there is the evidence of decreasing esp from their V°''^T, 
supporting this scenario. We do not wish to over-interpret 
the implications for ACDM at this point, since the num- 
ber of fainter galaxies in the sample is not yet statistically 
large, and most of them are not yet modelled in enough 
detail to exclude the potentially important effects of radial 
anisotropy at large radius. But this paper highlights the im- 
portance of detailed dynamical studies of a larger sample of 
galaxies, where the distribution of the dark matter as well 
as the stars may be used as clues to the formation histories 
of galaxies. 



ACKNOWLEDGMENTS 

We are very grateful to the referee. Dr. A. Bosma, for his 
comments which allowed us to improve the paper results. We 
thank James Bullock for providing his toy model code. NRN 
thanks M. Pannella for stimulating discussion on galaxy for- 
mation. NRN is receiving a grant from the EU within the 5th 
Framework Program (FP5) - Marie Curie Individual Fellow- 
ship. NRN, MC and MA acknowledge financial support by 
INAF-Project of National Interest. 



REFERENCES 

Alam S.M.K., Bullock J.S., Weinberg D.H., 2002, ApJ, 572, 
34 

Arnaboldi M. et al., 1996, ApJ, 472, 145 (A96) 
Arnaboldi M. et al., 1998, ApJ, 507, 759 (A98) 



16 Napolitano et al. 



Baes M., Dejonghe H., 2001, ApJ, 563, L19 

Bahcall N.A., Lubin L.M., Dorma V., 1995, ApJ, 447, L81 

Bender R., Burstein D., Faber S.M., 1992, ApJ, 399, 462 

Bender, R. & Nieto, J.L., 1990, A&A, 239, 97 

Benson A.J., Cole S., Frenk C.S., Baugh CM., Lacey C.G., 

2000, MNRAS, 311, 793 

Bernard! M. et al., 2003, AJ, 125, 1849 
Bertola F., Pizzella A., Persic M., Salucci P, 1993, ApJ, 
416, L45 

Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, 

J.R., 1986, ApJ, 301, 27 
Borriello A., Salucci P., Danese L., 2003, MNRAS, 341, 

1109 

Bosma A., van der Kruit P.C., 1979, A&A, 79, 281 
Bullock J.S., Kolatt T.S., Sigad Y., Somerville R.S., 
Kravtsov A.V., Klypin A. A., Primack J.R., Dekel A., 

2001, MNRAS, 321, 559 

Buote D.A., Jeltema T.E., Canizares C.R., Garmire G.P., 

2002, ApJ, 577, 183 

Busarello G, Capaccioli M., D'Onofrio M., Longo G., 
Richter G., Zaggia S., 1996, A&A, 314, 32 (B96) 

Capaccioli M., Caon N., D'Onofrio M., 1992, ESQ 
ESP/EIPC Workshop on Structure, Dynamics and Chem- 
ical Evolution of Early-type Galaxies (Elba), eds. J. 
Danziger, W. W. Zeilinger, and K. Kjar, ESO: Garching, 
43 

Capaccioh M., Cappellaro E., Held E. V., Vietri M., 1993, 

A&A, 274, 69 (C93) 
Capaccioh M., Napolitano N. R., Arnaboldi M., 2003, 

Sakharov Conference of Physics, Moskow, June 2002, 

[preprint astro-ph/0211323 
Chiosi C, Bressan A., Portinari L., Tantalo R., 1998, A&A, 

339, 355 

Cote, S., Carignan, C, & Sancisi, R., 1991, AJ, 102, 904 
Cote P., McLaughlin D.E., Cohen J.G., Blakeslee J. P., 

2003, ApJ, 591, 850 (C03) 

Danziger I.J., 1997, in Persic P., Salucci, P., eds, ASP Conf. 
Ser. Dark and Visible Matter in Galaxies, S. Francisco, 
ASP 

Davies, R. L. & Illingworth, G. D., 1986, ApJ, 302, 234 
de Blok, W.J.G., Bosma A., McCaugh S., 2003, MNRAS, 

340, 657 

Dekel A., Silk J. 1986, ApJ, 303, 39 

De Rijcke, S., Dejonghe, H., Zeilinger, W.W., & Hau, 

G.K.T. 2004, A&A, 426, 53 
de Vaucouleurs G., 1948, Ann. d'Astrophys., 11, 247 
de Zeeuw, P. T., et al., 2002, MNRAS, 329, 513 
Dirsch B. et al., 2004, AJ, 127, 2114 
Douglas N. et al., 2002, PASP, 114, 1234 
Dubinski J., 1998, ApJ, 502, 141 

Emsellem E., Dejonghe H., Bacon R., 1999, MNRAS, 303, 
495 (E99) 

Faber, S.M. et al., 1997, AJ, 114, 1771 
Forbes, D.A., Franx, M. & Illingworth, G.D., 1995, AJ, 109, 
1988 

Franx, M., van Gorkom, J. H., & de Zeeuw, T., 1994, ApJ, 
436, 642 

Fukugita M., Hogan C.J. & Peebles, P.J.E., 1998, ApJ, 503, 
518 

Gentile G., Salucci P., Klein U., Vergani D., Kalberla P., 

2004, MNRAS, 351, 903 

Gerhard O. et al., 2001, AJ, 121, 1936 



Gebhardt K. et al., 2003, ApJ, 583, 92 (G03) 

Gnedin O.Y., Kravtsov A.V., Klypin A.A., Nagai D., 2004, 

ApJ, submitted [preprint astro-ph/0406247 1 
Graham, A. W., Erwin, P., Trujillo, 1., & Asensio Ramos, 

A., 2003, AJ, 125, 2951 
Graham, A. W. & Guzman, R., 2003, AJ, 125, 2936 
Guzik J., Seljak U., 2002, MNRAS, 335, 311 
Hernquist L., 1990, ApJ, 356, 359 

Jimenez, R., Panter, B., Heavens, A., Verde, L., 2004, ac- 

cepetd on MNRAS [preprint astro-ph/0403294 
Jobin, M. & Carignan, C. 1990, AJ, 100, 648 
Kauffmann G. et al., 2003, MNRAS, 341, 33 
Keeton C.R., 2001, ApJ, 561, 46 
Koopmans, L.V.E. & Treu, T., 2002, ApJL, 568, L5 
Kormendy, J., 1977, ApJ, 218, 333 

Kronawitter, A., Saglia, R.P., Gerhard, O. & Bender, R. 

2000, A&As, 144, 53 (KOO) 
Kuntschner H., Lucey J.R., Smith R.J., Hudson M.J. & 

Davies R.L., 2001, MNRAS, 323, 615 
Lake, G., Schommer, R.A., & van Gorkom, J.H., 1990, AJ, 

99, 547 

Loewenstein, M. & White, R.E., 1999, ApJ, 518, 50 
Lokas, E.L., 2001, MNRAS, 327, L21 

Magorrian, J. & Ballantyne, D., 2001, MNRAS, 322, 702 
(MBOl) 

Mamon, G.A. & Lokas, E.L., 2004a, MNRAS, submitted 

[preprint astro- ph/0405466 
Mamon, G.A. & Lokas, E.L., 2004b, MNRAS, submitted 

[preprint astro-ph/0405491 
Marchesini D., D'Onghia E., Chincarini G., Firmani C, 

Conconi P., Molinari E., Zacchei A., 2002, ApJ, 575, 801 
Maraston, C, 1998, MNRAS, 300, 872 
Marinoni, C. & Hudson, M.J., 2002, ApJ, 569, 101 
Martimbeau, N., Carignan, C, & Roy, J.R., 1994, AJ, 107, 

543 

Mateo, M.L., 1998, ARA&A, 36, 435 

Matsushita, K., Belsole, E., Finoguenov, A. & Bohringer, 

H. 2002, A&A, 386, 77 
Mendez, R.H. et al., 2001, ApJ, 563, 135 (MOl) 
McGaugh S.S., Barker M.K., de Blok W.J.G., 2003, ApJ, 

584, 566 

Meza, A., Navarro, J.F., Steinmetz, M. & Eke, V.R., 2003, 

ApJ, 590, 619 
Mo, H.J. & Mao, S., 2004, MNRAS, 353, 829 
Moore B., Quinn T., Governato F., Stadel J., Lake G., 

1999, MNRAS, 310, 1147 
Morganti, R., Sadler, E.M., Oosterloo, T., PizzeUa, A., & 

Bertola, F., 1997, AJ, 113, 937 
Morganti, R., Oosterloo, T., & Tsvetanov, Z., 1998, AJ, 

115, 915 

Myers, S.T., et al., 1995, ApJL, 447, L5 

Napolitano N.R., Arnaboldi M., Capaccioh M., 2002, A&A, 

383, 791 (NAC02) 
Napolitano N.R., et al., 2004, lAUS, 220, 173 [preprint 

astro-ph/0310798 
Navarro J.F., Frenk C.S., White S.D., 1996, ApJ, 462, 563 

(NEW) 

Navarro J.F., Frenk C.S., White S.D., 1997, ApJ, 490, 493 
(NEW) 

Navarro J.F. et al., 2004, MNRAS, 349, 1039 
Nicto J.L., Bender R., 1989, A&A, 215, 266 



Mass-to-light ratio gradients in early-type galaxy haloes 17 



Oosterloo, T. A., Morganti, R., Sadler, E. M., Vergani, D., 

& Caldwell, N., 2002, AJ, 123, 729 
O'Sullivan, E. & Ponman, T.J., 2004a, MNRAS, 349, 535 
O'SuUivan, E. & Ponman, T.J., 2004b, MNRAS, 354, 935 
Padmanabhan N. et al., 2004, New Astronomy, 9, 329 
Pahre M.A., Djorgovski S.G., de Carvalho R.R., 1995, ApJ, 

453, L17 

Pahre M.A., Djorgovski S.G., de Carvalho R.R., 1998, AJ, 

116, 1591 
Pellegrini S., 1999, A&A, 351, 487 

Peng E.W., Ford H.C., Freeman K.C., 2004, ApJ, 602, 685 
(P04) 

Press W.H., Teukolsky S.A., Vetterling W.T. & Flannery 

B.P., 1992, Numerical Recipes in C, 2nd Ed. (Cambridge 

Univ. Press) 
Ravindranath S. et al., 2001, ApJ, 122, 653 
Richtler T. et al., 2004, AJ, 127, 2094 (R04) 
Romano, D., Silva, L., Matteucci, F., & Danese, L., 2002, 

MNRAS, 334, 444 
Romanowsky A. J., Kochanek C.S., 2001, ApJ, 553, 722 

(RKOl) 

Romanowsky A.J. et al., 2003, Science, 301, 1696 (R03) 
Rusin D., Kochanek C.S., Keeton C.R., 2003, ApJ, 595, 29 
Saglia R.P., Kronawitter A., Gerhard, O., Bender R., 2000, 

AJ, 119 153 (SOO) 
Saiz A., Dommguez-Tenreiro R., Serna A., 2004, ApJ, 601, 

L131 

Shen S., Mo H.J., White S.D.M., Blanton M.R., Kauffmann 
G., Voges W., Brinkmann J., Csabai I., 2003, MNRAS, 
343, 978 

Sommer-Larsen J., Gotz M., Portinari L., 2003, ApJ, 596, 
47 

Statler T.S., Dejonghe H., Smecker-Hane T., AJ, 117, 126 
(S99) 

Terlevich A.I., Forbes D.A., 2002, MNRAS, 330, 547 
Tonry J.L., Dressier A., Blakeslee J. P., Ajhar E.A., 

Fletcher A.B., Luppino G.A., Metzger M.R., Moore C.B., 

2001, ApJ, 546, 681 
Trager S.C., Faber S.M., Worthey G., Gonzalez J.J., 2000, 

AJ, 120, 165 

Tremblay B., Merritt D., Williams T.B., 1995, ApJ, 443, 
L5 (T95) 

Treu, T. & Koopmans, L.V.E., 2004, ApJ, 611, 739 
Trujillo I., Burkert A., Bell E. F., 2004, ApJ, 600, L39 
TuUy, R. B., 1988, Nearby Galaxies Catalog (Cambridge 
Univ. Press) 

van den Bosch F.C., Mo H.J., Yang X., 2003, MNRAS, 345, 
923 

van den Bosch E.G., Norberg P., Mo H.J. , Yang X., 2004, 

MNRAS, 352, 1302 
van der Marel R.P., 1991, MNRAS, 253, 710 (M91) 
Verolme E.K. et al., 2002, MNRAS, 335, 517 (V02) 
Walcher, C.J., Fried, J.W., Burkert, A., & Klessen, R.S., 

2003, A&A, 406, 847 
Wechsler R.H., Bullock J.S., Primack J.R., Kravtsov A.V., 

Dekel A., 2002, ApJ, 568, 52 
Weinberg, D.H., Dave, R., Katz, N., & Hernquist, L., 2004, 

ApJ, 601, 1 

Whitmore, B.C., Forbes, D.A., & Rubin, V.C., 1988, ApJ, 
333, 542 

Worthey, G., 1994, ApJs, 95, 107 

Wright, L.J., Ostriker, J. P., Naab, T. & Efstathiou, G., 



2004, submitted to MNRAS [preprint [astro-ph/OSlOSl^ 



